Constraints on the Brans-Dicke gravity theory with the Planck data 
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Based on the new cosmic microwave background (CMB) temperature data from the Planck satel- 
lite, the 9 year polarization data from the Wilkinson Microwave Anisotropy Probe (WMAP), and 
the baryon acoustic oscillation (BAO) distance ratio data from the Sloan Digital Sky Survey (SDSS) 
and 6 Degree field (6dF) surveys, we place a new constraint on the Brans-Dicke theory. We adopt 
a parametrization £ = ln(l + ~), where the general relativity (GR) limit corresponds to £ = 0. We 
find no evidence of deviation from general relativity. At 95% probability, —0.00246 < C < 0.00567, 
correspondingly, the region —407.0 < uj < 175.87 is excluded. If we restrict ourselves to the £ > 
(i.e. uj > 0) case, then the 95% probability interval is £ < 0.00549, corresponding to uj > 181.65. 
We can also translate this result to a constraint on the variation of gravitational constant, and find 
the variation rate today as G — — lAZt^'zr x 10 _13 yr _1 ( la error bar), the integrated change since 
the epoch of recombination is 5G/G = 0.0104lQ oJg7 (la error bar). These limits on the variation 
of gravitational constant are comparable with the precision of solar system experiments. 



I. INTRODUCTION 

The Jordan-Fierz-Brans-Dicke theory [lHH (hereafter 
the Brans-Dicke theory for simplicity) is the simplest ex- 
tended theory of gravity. In addition to the metric tensor, 
there is a scalar field in this theory, i.e. the Brans-Dicke 
field, which gives the effective gravitational constant 0- 
1 1 Oil - The action in the Jordan frame is 
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where is the action for the matter field. Here 4> is the 
Brans-Dicke field and uj is the Brans-Dicke parameter. In 
the limit of uj — > oo the Brans-Dicke theory is reduced to 
Einstein's general relativity theory. 

Solar system experiments have already put very strin- 
gent constraints on the Brans-Dicke mo del [Til. For 
example, the tracking data obtained from the Cassini 
mission gives uj > 40000 at the 2a level |l3|. Neverthe- 
less, it is still interesting to test the theory on cosmo- 
logical scales, especially because the Brans-Dicke theory 
may be regarded as an approximation for a number of 
scalar-tensor theories of gravity which have more signifi- 
cant effects on larger scales. The cosmic microwave back- 
ground (CMB) anisotropy can be calculated for a given 
theory, and the Brans-Dicke model may be tested with 
high precision (l4j . 

A number of limits on the uj parameter have been de- 
rived since WMAP released its data on CMB anisotropy. 
Acquaviva et al. report that uj > 80 at 99% level by 
combining the WMAP 1st year data and some ground 



or balloon based experiments and the large scale struc- 
ture data[lj|. Wu et al. [H, [TtJ excluded the region 
of -120 < uj < 97.8 by using the WMAP 5 year data 
and the Sloan Digital Sky Survey (SDSS) Luminous Red 
Galaxy (LRG) data. Considering only the possibility of 
uj > 0, Avilez and Skordis[l8j derived a limit of uj > 288 
at 95% confidence level, by combining the WMAP 7 year 
data and the data from the South Pole Telescope (SPT) 
and other small scale CMB experiments. When compar- 
ing these different results, one should note that the limits 
obtained depend very much on the parameterization and 
prior used, see the next section for discussion. 

The precision of cosmological observations are being 
improved steadily. The WMAP grou p have published 
the data of 9 years of observation pjj [2(| , and recently, 
the Planck collaboration published their observational re- 
sults [2l[. In addition to the CMB observations, there are 
also much progress in redshift surveys of galaxies. Recent 
surveys such as the SDSS-III Baryon Oscillation Spec- 
troscopic Survey (BOSS) 1 and 6dF 2 have measured the 
power spectrum of the large scale structure at different 
redshifts, and obtained cosmic distances from the baryon 
acoustic oscillation features. 

In this paper, we update the constraint on the Brans- 
Dicke model by using the new CMB data 3 , including the 
Planck temperature anisotropy 22( and the the WMAP9 
CMB polarization data[lj|. Following the Planck collab- 
oration [23[, in addition to the CMB data, we also use 
the BAO data from the SDSS 0-11] and 6dF H3] galaxy 
redshift surveys. 
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II. METHODS 

For convenience, we introduce a dimensionless field 
tp = G<p where G is the Newtonian gravitational con- 
stant, then the Einstein equations are 
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where T)i™^ is the stress tensor for matter. The equation 
of motion for ^ is 
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For G to be consistent with the Cavendish experiment, 
The value of ip at present day should be 

2w + 4 
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We follow the calculation method described in (16j . 
in which we developed the covariant and gauge-invariant 
formalism of cosmological perturbation theory in the case 
of Brans-Dicke gravity, and apply the method to calcu- 
late the angular power spectra of CMB temperature and 
polarizations, as well as the power spectrum of large scale 
structure (LSS). 

Given a cosmological model, the angular power spec- 
tra of CMB temperature and polarization and the matter 
power spectrum can be calculated, for example with the 
publicly available code CAMB [28| . In order to constrain 
the cosmological parameters with the observational data, 
we use the publicly available CosmoMC code [2j|, which 
uses the Markov Chain Monte Carlo (MCMC) method to 
explore the parameter space, with a modified CAMB code 
developed by Wu et al. in Ref JT3]- We use the latest 
CMB data published by the Planck team [22|]. Accord- 
ing to our previous analysis fl6j . the small scale (high-?) 
anisotropy is affected more by the Brans-Dickc gravity, so 
the more precise measurements of Planck on small scales 
(up to I ~ 2500) should help greatly. For the low-Zs, we 
also include the TE and BB power spectrum estimated 
from the polarization map of WMAP9, though the latter 
does not provide much distinguishing power at present. 

We also combine the BAO data from large scale struc- 
ture surveys, including the SDSS DR7 |j, 1|, BOSS 
DR9 [Hi and 6dF [2j. The BAO surveys measure the 
distance ratio 
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where r s (zdrag) is the comoving sound horizon when 
baryons became dynamically decoupled from photons 
(the baryon drag epoch) and Dy(z) is the combination 
of angular-diameter distance, Da{z), and the Hubble pa- 
rameter, H(z) 



D v {z) = [{l + zfD\(z) 



H(zY 



11/3 



(0) 



We follow the choice of BAO data set "SDSS DR7 + 
BOSS DR9 + 6dF" in the Planck analysis 0. This 
includes two of the most accurate BAO measurements, 
and minimizes the correlations between the galaxy sur- 
veys, as the two surveys have widely separated effective 
redshifts. 

The derived limits depend on which parameterization 
is selected and how the priors are set. As the experiments 
so far all favors the GR case and the Brans-Dicke param- 
eter is stringently constrained, it is more convenient to 
take the GR as the null case, and have a parametrization 
in which the tested parameter vanishes for GR. In prac- 
tice a flat prior on a finite range is usually assigned to 
the parameter. For example, Ref.[l5[ and Ref.[l8[ con- 
sidered flat prior on ln[l/w] (though the exact parameter 
they used differs slightly). A limitation of this choice is 
that it could not treat the to < case. Ref. [l8| argued 
that if tu < — 3/2 the Brans-Dickc field would be a ghost 
field, and they will therefore consider only positive u>. 
However, we would rather err on the conservative side, 
and use a more general form of parameterization which 
allows negative u. Indeed, at present there are many 
phantom dark energy 



also ghost-like. In Ref. 



301 models in which the field are 
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C = ln(l + -) , 



(7) 



in the present paper we will also adopt this parameter- 
ization. This parameter has the nice property that as 
( — >• 0, the Brans-Dicke theory reduces to Einstein grav- 
ity, and it is easy to obtain limits on both the negative 
and positive value of ui. We choose the same initial range 
[-0.014,0.039] as Wu ct al. in 0, which is convenient 
for computation, while at the same time the final con- 
straint is not very sensitive to this range, since at the 
edge of the prior range the likelihood is very small. In 
fact, if one wishes to consider only positive values of ui, 
we can also do that by simply restricting the range of the 
prior to < C < 0.039. 

As pointed out by Ref. (l8| . comparing with their prior, 
our flat prior on £ penalizes large ui, and hence for the 
same data set a "weaker" limit on ui would be obtained 
for our choice. We do not see a good theoretical reason 
to favor one prior on ui over the other, but our choice 
is again in agreement with our general philosophy of be- 
ing conservative on constraining models. We remind the 
reader to notice the effect of the prior when comparing 
results obtained in different papers. 

We also obtain limits on the following basic or de- 
rived cosmological parameters: ftbh 2 , £l c h 2 , 0, r, 
n s , ln(10 10 Aj), Age/Gyr, erg, z le and Hq. Here J7a is the 
dark energy density today, fib is the baryon density to- 
day. f2 c is the cold dark matter density today. 6 is the 
angular scalar of the sound horizon at last-scattering, 
r is the Thomson scattering optical depth due to the 
rcionization. n s is the scalar spectrum power-law index. 
In(10 10 j4 s ) is log power of the primordial curvature per- 
turbations. Age/Gyr is the age of the universe, eg is 
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the rms matter fluctuations today in linear theory. z IC is 
the redshift at which universe is half reionized. H is the 
Hubble constant. 



III. RESULTS 

FiglD shows the 1-D marginalized distribution for the 
Brans-Dicke parameter £. The curve which is labelled as 
"Planck + WP" shows the CMB-only result, for which 
the temperature data from Planck and the polarization 
data from WMAP9 is used. The curve which is labelled 
as "Planck + WP + BAO" combined the CMB data with 
the BAO observation data from SDSS DRY, BOSS DR9 
and 6dF. For comparison, we also plot in this figure the 
result obtained in our previous work [l7j , which used the 
WMAP5 data and the matter power spectrum from the 
SDSS LRG survey. 




C 

FIG. 1: The one dimensional likelihood distribution for £. 
"Planck + WP" denotes the result of using Planck temper- 
ature data as well as WMAP9 polarization data. "Planck 
+ WP + BAO" denotes the combined constraint with BAO 
data 03 HH [2(| [27|]. We also plot the result in previous 
work "WMAP5 + SDSS LRG", which is using CMB 

temperature and polarization data from WMAP5, combined 
with matter power spectrum measured with the luminous red 
galaxy (LRG) survey of the SDSS. 



Comparing with the previous work, especially Wu et 
al.[l7|. which used the same parameterization, the con- 
straints become stronger, and the new data favors a 
slightly more positive value of £. The likelihood of the 
CMB-only data looks quite Gaussian. Because of the 
high angular resolution of Planck data, the CMB-only 
data can already give better constraints than before. 
With the BAO data, the constraint is further tightened. 
The BAO distances estimated from galaxy surveys play 
the same role as the matter power spectrum in distin- 
guishing the different models. 

For the CMB-only case, we find the 68% and 95% in- 



tervals arc 

- 0.247 x 1CT 2 < C < 1-080 x 1CT 2 (68%) ; (8) 
-0.855 x 10~ 2 < C < 1-716 x 10~ 2 (95%) . (9) 

These correspond to 

uj < -405.36 or uj > 92.09 (68%) ; (10) 
uj < -117.46 or uj > 57.78 (95%) . (11) 

We see that with the CMB data alone, the constraint is 
still relatively loose. 

Addition of the BAO data helps to break the parameter 
degeneracy, and much stronger limits are obtained. For 
the CMB+BAO case, we find the 68% and 95% bounds 
are 

- 0.046 x 10~ 2 < C < 0.366 x 10" 2 (68%) ; (12) 
-0.246 x 10~ 2 < C < 0.567 x 10~ 2 (95%) , (13) 

which correspond to 

uj < -2174.41 or uj > 272.72 (68%) ; (14) 
uj < -407.00 or uj > 175.87 (95%) . (15) 

If we restrict ourselves to the case of £ > 0, or equiv- 
alently uj > 0, then for the CMB only case the 68% and 
95% bounds are 

< C < 0.895 x 10~ 2 (68%) ; (16) 
< C < 1-645 x 10~ 2 (95%) , (17) 

corresponding to 

uj > 111.23 (68%); uj > 60.29 (95%) . (18) 

For the CMB+BAO case, 

< C < 0.296 x 10~ 2 (68%) ; (19) 
< C < 0.549 x 10" 2 (95%) , (20) 

corresponding to 

uj > 337.34 (68%); uj > 181.65 (95%) . (21) 

We see that the constraints are only slightly different 
from their respective positive bounds where £ < 0(uj < 0) 
are allowed, even though the a prior allowed parameter 
space is smaller. This shows that the negative C solutions 
fit about also very well, so reduction of parameter space 
does not significantly improve the constraint. 

As discussed in the last section, these limits depend 
on the parameterization and prior adopted. The results 
presented here applies to the parameter £, even though 
we also quoted limits on uj since that's what appeared in 
the Brans-Dicke theory. This parameterization is more 
"conservative" , so our limits appeared to be "weaker" 
than Ref.[l8| even though we have used the newer and 
more precise Planck data. 
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FIG. 2: The two dimensional contour for £ against J1a- 



Fig [2] shows the two dimensional contours for £ against 
f^A • If only the CMB data from Planck is used, the con- 
straint already become stronger, but there is significant 
degeneration between £ and J7a- The BAO data can help 
to break the degeneration and give much stronger con- 
straints. The center of the contours shifted somewhat 
from the center of our previous results, this is the same 
trend as seen in the fitting of the standard cosmological 
model, for is lowered, but we see that the shift on the 
center of £ is small. 

Next we examine how the constraints on other cosmo- 
logical parameters are affected if we consider the Brans- 
Dicke gravity. FigJ3]shows the one dimensional likelihood 
for some cosmological parameters. In this plot, we show 
the result with £ fixed to 0, labelled as "GR", in this 
case the Brans-Dicke theory is reduced to the standard 
ACDM model with Einstein's General Relativity. It is 
obvious from the figure that when the BAO data are com- 
bined, the constraint is much tighter than the case with 
CMB data only. However, for most parameters, the like- 
lihood distribution of the GR case and the Brans-Dicke 
case is very similar, the shifts in the best fit parameters 
(peak value of the likelihood) are small, and the differ- 
ences in the width of the likelihood are also relatively 
small, showing that the addition of the Brans-Dicke pa- 
rameter does not significantly affect the uncertainty in 
other cosmological parameters. The most affected basic 
parameters are Ho, £l c h 2 , and as, while for the derived 
parameters the uncertainty on the cosmic age is much 
larger. 

We also plot the two dimensional contours in FigJS] 
The Planck data can give accurate measure on tl^h 2 
without other additional data, thanks to the well mea- 
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FIG. 3: The one dimensional likelihood for cosmological pa- 
rameters. The red lines with label "Planck" represent the re- 
sult with Planck temperature data and WMAP9 polarization 
data; the blue lines with label "Planck + BAO" denote the 
result combined with BAO data; the green line labelled "GR" 
is the result of fixing £ = 0, in which case the Brans-Dicke 
gravity reduces to Einstein theory. 



sured peaks in angular power spectra. The degeneracy 
between £ and Qbh 2 , as well as £l c h 2 are quite limited and 
the uncertainty of fitting is reduced. The degeneracy be- 
tween £ and f2 m , £l\ are broken by BAO data, which 
could also be seen in one dimensional likelihood distribu- 
tion, f2 ro and Q\ change greatly after introducing BAO 
data. 

The best-fit values and 68% marginalized error are 
shown in Table [TTJ For comparison, in Tabic |TT] we also 
list the results for Einstein gravity given in Ref . [22| . The 
Einstein result is constrained by using Planck low-/ and 
high-/ data, as well as WMAP9 polarization data and 
BAO data, which is the same as our data set. Compar- 
ing with the WMAP data, the Planck data favors lower 
Oa and lower Hq, in the standard ACDM model fitting. 
This trend was noted by the Planck team and also found 
in our model. 

There are some slight differences between our results 
with £ fixed to and the result published by the Planck 
team [23}. The differences come mainly from different 
setting of parameters in CosmoMC. In our fitting, in or- 
der to focus on the Brans-Dicke parameters, we ignored 
the effect of massive neutrinos, and fixed the neutrino 
number. The best-fit values of cosmological parameters 
in our Brans-Dicke model are consistent with ACDM 
model in Einstein theory. 

We can also test the variation of the gravitational 
constant G in the context of Brans-Dicke theory. We 
added two derived parameters in the MCMC code, i.e. 
G/G = —<p/<p, which is the change rate of gravita- 
tional constant at present, and SG/G = (G rec — Gq)/Gq, 
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which is the integrated change of gravitational constant 
since the epoch of recombination. The one dimensional 
marginalized likelihood is shown in Fig 4(a) and Fig |4(b)| 
The likelihood functions are fairly close to the Gaussian 
form. We can take the 68% limit as corresponding to the 
la error for these measurements. 

For the CMB only case, the best-fit values are 

G/G = -0.4617 x 10~ 12 , SG/G = 0.0318 

and the 68% marginalized limits are 

- 1.1970 x 10~ 12 < G/G < 0.4597 x 10" 12 ; (22) 

-0.0197 < SG/G < 0.0835 (23) 

For the CMB+BAO case, 

G/G = -0.1417 x 10 -12 , 5G/G = 0.0104 
and the 68% marginalized limits are 

- 0.4082 x 10~ 12 < G/G < 0.0663 x 10~ 12 (24) 

-0.0037 < SG/G < 0.0290 (25) 

We list the constraints on G/G with different methods 
in Table HI Though model-dependent, our cosmological 
constraints are now comparable in precision with other 
methods, including the solar system experiments. 



TABLE I: Constraints on the rate of variations of gravita- 
tional constant. The errors are la unless otherwise noted. 



larization data [19| to constrain the Brans-Dickc theory. 
In addition to the Planck data, we also use the BAO 
data from the SDSS DR7 0, BOSS DR9 0I and 



G/G [HT^yr- 1 ] 


Method 


2±7 


lunar laser ranging [31] 


0±4 


big bang nucleosynthesis [321 [33] 


0± 16 


helioseismology [341 


-6 ±20 


neutron star mass [35] 


20 ±40 


Viking lander ranging [361 


40 ±50 


binary pulsar [37] 


-96 ~ 81 (2a) 


CMB (WMAP3) [38] 


-17.5 ~ 10.5 (2a) 


WMAP5±SDSS LRG [17J 


-1.42121? 


Planck±WP±BAO (This paper) 



6dF [27J , which help to break parameter degeneracy. 



We use the parameterization £ = ln(l + —) introduced 
in Ref.[l7[, for which the GR limit is achieved when 
C — > 0, (\u\ — > oo). This parameterization may be more 
"conservative" than some of the other parameterizations, 
so the limit we derive may also appear "weaker" than 
given in some of the other works. The readers should 
note this when comparing the results given in different 
works. 

We obtained constraints by using the CMB data (re- 
ferred to as CMB-only). The 68% and 95% bounds are 
given Eqs.([8|)- (|lip . By combining the BAO data, we 
obtain stricter constraints, which are given in Eqs. ([T2")) - 
(TT5l) . We also considered the bounds obtained if ( > 
(or equivalently uj > 0) is assumed to be positive, these 
are given in Eqs. ([T6"|) - (f2"Tj) . We do not detect any sig- 
nificant deviation from Einstein's general theory of rel- 
ativity, and the constraint on the Brans-Dicke model is 
tightened compared with previous results. 

We examined the distribution of other cosmological pa- 
rameters. For most parameters, the best fit values and 
measurement errors are not altered much by the intro- 
duction of the Brans-Dicke gravity. The most affected 
parameters are Hq, f2 c /i 2 , and as, and the derived pa- 
rameter of cosmic age. 

Finally, the variation of the gravitational constant in 
the Brans-Dicke model are also constrained, the results 
are given in Eqs. (f2"2"]) - (|2"5|) . and also summarized in Table 
HI These constraints are model-dependent, nonetheless, 
it is remarkable that the limits obtained are comparable 
with the constraints from the highly precise solar system 
experiments. 
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TABLE II: Summary of cosmological parameters and the corresponding 68% intervals. The "Planck + WP" column lists the 
result of using temperature map from Planck and polarization map from WMAP9; The "Planck + WP + BAO" column lists 
the result with BAO data combined; We also list the result using the same data as "Planck + WP + BAO", but fix £ = in 
the "Planck + WP + BAO with £ = 0" column, that Brans-Dicke reduces to Einstein theory. The last column is the result 
form Planck team in Ref . [22l| . 
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FIG. 5: The two dimensional contour for cosmological parameters. 



